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ON EXTREMA OF STABLE PROCESSES 

By Alexey Kuznetsov^ 

York University 

We study the Wiener-Hopf factorization and the distribution of 
extrema for general stable processes. By connecting the Wiener-Hopf 
factors with a certain elliptic-like function we are able to obtain many 
explicit and general results, such as infinite series representations and 
asymptotic expansions for the density of supremum, explicit expres- 
sions for the Wiener-Hopf factors and the Mellin transform of the 
supremum, quasi-periodicity and functional identities for these func- 
tions, finite product representations in some special cases and iden- 
tities in distribution satisfied by the supremum functional. 

1. Introduction. Stable processes are probably the most fascinating mem- 
bers of the vast family of Levy processes. They enjoy the scaling prop- 
erty, which states that the process {X^ : i > 0} has the same distribution as 
{c^/'^Xt :t>0}. Another well-known process which also enjoys this property 
is Brownian motion, for which we can compute explicitly the distribution 
(in many cases even the joint distribution) of many important functionals; 
an extensive collection of these facts can be found in [8] . In the case of sta- 
ble processes the situation is much less satisfactory: we do not have general 
explicit results even for such a basic functional as the supremum. However, 
there does exist an "almost" general result: Doney [11] has obtained closed- 
form expressions for the Wiener-Hopf factors for a dense set of parameters. 
These results have inspired the main idea for our current work. 

In [11] it was proven that if parameters a and p = P(Xi > 0) satisfy 

(1.1) p + k = - 

a 

for some integers k and I, then there exists a fully explicit representation for 
the Wiener-Hopf factors, given as a ratio of two finite products of length 
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\k\ and |/|. It is easy to see that when a is irrational, relation (1.1) defines 
a countable set of values of p, which is dense in [0, 1]. Thus, assuming that 
a is fixed and irrational, we have a function with explicit representation for 
a dense countable set of values of the parameter p. This situation reminds 
one of an analogous property of elliptic functions: when a certain parameter 
is a rational number there often exists an explicit expression, or at least 
a reduction to some simpler functions. Therefore, given results obtained in 
[11], it seems highly probable that there should exist a connection between 
extrema of stable processes and elliptic functions. The goal of this paper is 
to exhibit this connection and as a consequence to obtain a remarkably large 
number of explicit results on extrema of stable processes. 

First we will present some definitions and notations. We will always 
use the principal branch of the logarithm, which is defined in the domain 
|arg(z)| < TT by requiring that ln(l) = 0. Similarly, the power function will be 
defined as = exp(aln(z)) in the domain |arg(2:)| < vr. We will work with a 
stable process Xt, which is a Levy process with the characteristic exponent 
^'(z) = — ln(E[exp(z2;Xi)]) given by 

(1.2) ^{z) = c\z\'^(^l-il3tan(^^^sgn{z)^, z£R, 

where c> 0, a G (0,1) U (1,2) and /3 G [-1,1] (see [5] and [11]). The fact 
that the characteristic exponent is essentially a power function of order a 
implies the scaling property, which states that processes {Xd - t >0} and 
{c^/°Xj :t>0} have the same distribution. 

Following [11] (and [26]) we rewrite (1.2) in a more convenient form. First 
of all, c is just a scaling parameter, thus without loss of generality we can 
assume that = 1 + /3^tan(^)^. Next we introduce a parameter 

2 1 / ^ / vra 

(1.3) 7 = -tan^M -/3tan( — 

and rewrite (1.2) as 

(1.4) ^{z) = e-^/2|z|n|,>o} + e--^/2|z|n|,<o}- 

Another important parameter is p = P(Xi > 0), which was computed in 
closed form in [26] 

(1.5) p=l(l-2)=i + -L.a„-'(,8ta„(f)). 

One can see that parameters (a, /?) belong to the set 

{a G (0, l),/3 G (-1, 1)} U {a = l,/3 = 0} U {a G (1,2), /3 G [-1, 1]}, 
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if and only if parameters {a,p) belong to the following set: 
A = {ae{0,l),pe (0, 1)} u {a = 1, p = i} 

(1.6) 

U{aG (l,2),pG [1- a~^,a~^]}. 

We will call A the admissible set of parameters and will parametrize the 
stable process Xf by {a,p). Note that we exclude the case when a S (0, 1) 
and p = 1 {p = 0}, as in this case the process Xt {—Xt} is a subordinator 
and the Wiener-Hopf factorization is trivial. In the limiting case as a — )• 2~ 
and p — 7- 1/2 we obtain Brownian motion Xt = \plWt- When a € (1, 2) and 
p = 1 — [p = a~^) the process Xt is spectrally positive (negative). In 
this case we have complete information about the Wiener-Hopf factorization 
and the distribution of extrema due to the work of Bingham [7], Doney [12], 
Bernyk, Dalang and Peskir [4] and Patie [23]. 

The following family of stable processes, which was introduced in [11], 
will be very important in our paper. 

Definition 1. For k,l ^Z, define C^^i as the class of stable processes Xt 
with parameters (a,p) G A satisfying relation (1.1). 

Note that if {a,p) satisfy relation (1.1), then a G Q if and only if p G Q. 
When a ^ Q there exists a unique pair of integer numbers /c, I such that 
relation (1.1) holds. If a = ^ for some coprime integers m,n and Xt G Ck^u 
then relation (1.1) holds for any pair {k,l) = {k + jn,l + jm), j G Z. In this 
case we will assume that <k <n and 1 <l <m. Another important fact 
is that the process Xt with parameters {a,p) belongs to C^^i [and a G (1,2)] 
if and only if the process Xt with parameters {a,p) = {ap, l/a) belongs to 
C-i-k- This property, which can be easily verified using relation (1.1), will 
appear in many formulas throughout this paper and will be very useful for 
us later. Another important property is that Xt belongs to Ck^i if and only if 
the dual process Xt = —Xt belongs to C^k-i-i- This can be easily checked 
using (1.1), as the dual process has parameters (a, 1 — p). 

Next, we define the supremum and infimum processes 

St = sup{Xu ■■0<u<t}, It = mf{Xu ■.0<u<t}. 

We introduce a random variable e{q) ~ Exp(g) (exponentially distributed 
with parameter q > 0) which is independent of the process Xt- We will use 
the following standard notation for the characteristic functions of 5'e(|j) and 
/e(g), also known as the Wiener-Hopf factors: 

(/.+ (z) =E[e*^^-W], </.-(z) =E[e*^^<5)], z G M. 
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Note that since the random variable S(,(q) {/e(g)} is positive {negative}, 
function (p'^iz) {(pqiz)} can be extended analytically into domain Im(2:) > 
{lm{z) <0}. 

To reduce the complexity of the problem, we note first that it is enough to 
study the positive Wiener-Hopf factor (j)^{z), as the corresponding informa- 
tion about 4>~{z) can be obtained by considering the dual process. Second, 
the scaling property implies that 

0+(z) = </.+ (zg-V"); 

therefore we only need to consider the case (7 = 1. This justifies our choice 
of the following function as the main object of study in this paper: 

(1.7) (/>(z) = (/>(z;a,p) =E[e-'^'^'=(i)], Re(z) > 0. 

This paper is organized as follows: in Section 2 we establish a connection 
between (i){z) and a certain elliptic- like function F{z;t), and we study var- 
ious analytical properties of the latter. In Section 3 we derive an explicit 
expression (in terms of the Clausen function) for the function (f>{z) when 
a is rational. In Section 4 we study the Wiener-Hopf factorization in the 
general case: we express the Wiener-Hopf factor (f>{z) in terms of the q- 
Pochhammer symbol for Im(a) > 0, rederive some of the results obtained 
in [11], express (j){z) in terms of the double gamma function and give a 
complete description of the analytical properties of (j){z). In Section 5 we 
obtain a series representation for ln{(j){z)) with very interesting convergence 
properties and establish several functional equations satisfied by (piz). In 
Section 6 we study the Mellin transform of the supremum Si: we obtain 
two functional equations and establish quasi-periodicity of this function, ex- 
press it in terms of the double gamma function and as a finite product when 
S Ck,i and present two distributional identities satisfied by the supremum 
functional. Finally, in Section 7 we derive several convergent and asymptotic 
series representations for the probability density function of the supremum 
Si. 

2. Connecting the Wiener Hopf factors with elliptic functions. Our main 
tool in studying the Wiener-Hopf factor (j){z) will be a certain function 
F(z;t), which has many properties similar to elliptic functions, but first we 
need to define two types of domains in the complex plane. 

Definition 2. Assume Im(r) > 0. Define 

S{t) = {z G C : |Re(zf )| < ^lm(r)}, 
V{t) = S{t) r\{z€C: \lm{z)\ < 7rlm(r)}. 
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One can see that 5(r) is a strip in the complex plane which contains 
z = 0, and V{t) is a parallelogram if Re(r) 7^ 0. In the case when Re(r) = 
we have 5(r) = {z: |Im(2:)| < vr}, therefore 5(r) becomes a horizontal strip 
which does not depend on r, while parallelogram V{t) degenerates into the 
horizontal strip {z G C : |Im(2;)| < 7rmin(Im(r), 1)}. Domain V{t) satisfies 
the following important property: z G V{t) if and only if (—iz/r) G V{—1/t). 
Also note that 5(ar) = S{t) for a G M"*", therefore 5(r) depends only on the 
argument of r. 

Definition 3. Assume Im(r) > and z G 5(r). Define 

(2.1) F{z;r) ^ 



(l + e^+*™')(l + e^)' 

The integrand in the definition (2.1) converges to zero exponentially as 
X — )• +00, condition Im(T) > ensures that the same it true as x — t- —00 and 
condition z G S{t) guarantees that exp(2; + irx) 7^ —1 for 2; G M. Therefore 
the integral in (2.1) converges absolutely, and the function F{z;t) is well 
defined. It is also clear that this function is analytic in both variables {z, r) G 
as long as Im(r) > and z G S{t). 

Proposition 1. Assume {a,p) G A and |arg(z)| < 7rmin((l - p),^). 
Then 



(2.2) 



-^ln(0(z)) = J-[F{7Tip-ln{zy,ia-'] 
az Zttiz 



F{—7rip — ln{z);ia ^)]. 

Proof. We start with the following general integral representation for 
(/)+(z), which was first derived in [22] (Lemma 4.2) (another proof was given 
in [19], Theorem 1(b)): 

(2.3) ln((/.+ (z)) = -^ /inf ^) ^ , Im(z) > 0. 

As was shown in [22] , this formula is valid for any Levy process Xt , provided 
that the integral 

du 

converges for some e > 0, which is true in the case of stable processes due to 

(1.4) . Integral representation (2.3) is in fact equivalent to Darling's integral 
(see [10, 16]) which was the main tool used in [11]. 
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Next, we use our definition (j){z) = (l)f{iz) and formulas (1.4) and (2.3) to 
find that 



(2.4) 



ln(0(.)) = -- 



ln(l + e'^'T/^ 



u 



■ du 







+ 



u[u — IZ) 
ln(l + e-'^*'^/2^°) 







U[U + IZ 

We assume first that z G M'^ and obtain 

ln(l +e''^^/2^"^ 



du 



Reiz) > 0. 



d_ 

dz 







(2.5) 



u{u — iz) 

ln(l + e'^*T/2^") 
{u — izY 



du 



du 



d 



ln(l + e"^^/2^°) ^ \{u - izY^\ du 



-a 



,Q-1 



where in the last step we have applied integration by parts. Changing the 
variable of integration u = e~^/° in the last integral in (2.5) we have 



(2.6) 



a 



u 



a-l 



(e-'^*'>/2 + ii°)(M-iz) 



du ■■ 



dy 



[ervl^ - iz){\ + e-'^»T/2+y) ' 



The next step is to shift the contour of integration M 1— )■ M + ^ in the inte- 
gral in the right-hand side of (2.6). This step is justified, since the integrand 
decays exponentially as |Re(y)| — )■ 00, and the inequality I7I < a [which fol- 
lows from (1.3)] guarantees that the integrand is analytic in the horizontal 
strip Im(y) < |^|. Thus, shifting the contour of integration and performing 
one final change of variables y = x -|- ^ we obtain 

dy 



{e~yl°' - iz){l + e-'^n/2+j/) 

dy 



(2.7) 



h^i7/2 {e-y/'^ - iz){l + e--^T/2+?/) 
dx 



if dx 

z 



'X(i 



_|_ g7ri/2(l-7/Q)-ln(z)-x/o')^2 -|- e^) 
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Combining (2.5), (2.6), (2.7) with the definitions of F(z;t) and p we con- 
clude that 



In the case z G equation (2.2) follows by taking the real part of both sides 
of the above equation and using (2.4). We can extend this result by analytic 
continuation into domain |arg(z)| <7rmin((l — p),^), since the right-hand 
side of (2.2) is analytic in the region |arg(2;)| < 7r(l — p) and due to (2.4) 
ln((j){z)) is well defined and analytic in the domain |arg(z)| < f ■ □ 

At this point we would like to make several remarks which will be impor- 
tant later. First of all, the integral representation (2.4) allows us to define the 
function (p{z;a,p) for all positive a, even though for a > 2 these functions 
do not have an immediate probabilistic interpretation. Second, and most im- 
portantly. Proposition 1 allows us to consider (p{z;a,p) as a function which 
is analytic in all three variables {z, a, p). We do not need to specify explicitly 
the domain in C'^ where this function is analytic; for our purposes it is enough 
to know that for e > small enough the function <f)(z;a,p) is analytic in all 
three variables in the domain {|2 — 1| < e, |arg(a)| < e, |p— || < e} C [this 
follows immediately from (2.2) and the domain of analyticity of F{z;t)]. 

In the next theorem we collect several properties of the function F{z;t). 
As we will see later. Proposition 1 will allow us to translate each of these 
properties into an important statement about the Wiener-Hopf factor ^(z). 

Theorem 1. Assume Im(r) > 0. 




(i) For zeVir) 



(2.8) 




(ii) For z £ S{t) and n>2 



(2.9) 




(iii) ForzeV{T) and < e < min(f, Im(-^)) 



(2.10) 




(iv) For zeS{T)n-S{T) 



(2.11) 



F{z;t)=F{-z;t)--. 
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(v) Assume Re(T) / and define 6 = sign(Re(r)). For z £ 5(r) 

27ri 



(2.12) 



A;>0 



l+exp(z + 27r(^(A; + l/2)r) 
27rr-i 



+ 



1 + exp(iz/r + 2TT6/T{k + 1/2)) 
(vi) Assume Re(T) / 0. For z £ S{t) n {z : Re(z) > 0} 



(2.13) 



F(z;r) 



fc>i 



-kz 



-ikz/r 



sinh(7rA:T) r sinh(7rA:/T) 



Proof of Theorem 1. (i) Assume z > and Re(r) = 0. Then ir~^ G 
M^, and performing the change of variables x = —^y + ^ in the integral 
defining F{z;t) [see (2.1)] we immediately obtain identity (2.8). The general 
case can be obtained by analytic continuation, as both sides of (2.8) are 
analytic in the domain {Im(r) > 0,z G V{t)}. 

(ii) We start with the following identity: 



(2.14) 



1 1 

1 + n ^ 1 

fc=0 



_j_ y^g7ri/n(n— 2fc— 1) ' 



which is just a partial fraction decomposition of the rational function (1 + 
. Again, let us assume that z G M and Re(r) = 0. Applying this identity 



n\-l 



w 

to (2.1) we obtain 

dx 



(1 + e^+*"^^)(l + e^) 



1 """^ f 
n ^ X (I 



dx 



_j_ ^z/n+ni/n{n—2k—l)+iTx'j^'^ _j_ gX^ ' 



which is equivalent to (2.9). Using property 5(nr) = S{t) and the fact that 
z G 5(r) implies ^(-z i Tri{n — 1)) G 5(t), we find that both sides of (2.9) 
are analytic in {Im(r) > 0, z G 5(r)}, thus statement (ii) follows by analytic 
continuation. 

(iii) Assume z G M and Re(r) = 0, the general case can be established as 
usual by analytic continuation. We rewrite (2.1) as 



-(l+iT)x/2 



(2.15) F(z;r) = -e"^/2 / — , , , , — ---dx. 

^ ^ ^ ' ^ 4 cosh(x/2)cosh((z + irx)/2) 

Using formula 3.511.4 from [17] we find that for y G M and |t| sufficiently 
small 



-{l+iT)x/2+ixy 

cosh(x/2) 



■ dx 



sinh(7r(y — r/2)) ' 



-ixy 
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Applying Plancherel's theorem to (2.15) and using the above Fourier trans- 
form identities, we find that 

F(z;t) = / dy. 

2r cosh(7riy/r) sinh(7r(y - t/2)) 

To finish the proof, we shift the contour of integration in the above integral 
M I— )• M + I — ze (where e > is a sufficiently small number) and perform the 
change of variables y = + ^. 

(iv) Let f{x) be the integrand in (2.10). In order to derive identity (2.11) 
we start with the integral representation (2.10), shift the contour of inte- 
gration M -|- ie I— 7- M — ie taking into account the residue at x = and finally 
change the variable of integration x = —y 

2F{z;t)= [ fix)dx = -27riRes,=of{x)+ [ f{-y)dy 

= -27ri Res^=o f{x) + 2F{-z; r). 



z 

TVT ' 



The residue at zero of f{x) is easily seen to be 

(v) Let us fix n > 1 and assume Re(r) > 0. By shifting the contour of 
integration in (2.1) M i— t- M — 27rm (while taking care of the residues) and 
then changing variable of integration x = y — 27rin, we find 



F(z-T)=y 

^ l-hexp(z + 27r(fc + l/2)r) 

27rT~i 



^^^^ - + exp(iz/r + 27r/r(A; + 1/2)) 
dy 



~^ ' (1 + e^+2"'^^+«^2/)(l + e?') ' 

where n* = max{fc:Im(— ^ — ^^{k + ^)) < 27rn}. As n — )• +oo we have 
|g2n7rr| _^ which implies that the integral in the above equation converges 
to zero, and since n* — )• oo we obtain series expansion (2.12). If Re(r) < the 
proof is identical, except that now we need to shift the contour of integration 
into the upper half-plane: M i— t- M + 27rm. 

(vi) The proof is essentially the same as in part (v): we start with the 
integral representation (2.10) and shift the contour of integration into the 
upper-half plane while taking care of the residues. The details are left to the 
reader. □ 



3. Finite product representation when o: is rationaL As a first illustra- 
tion of the power of Proposition 1 and Theorem 1 we will derive an explicit 
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expression for (j){z) when a is rational and p is completely arbitrary. This 
expression involves the Clausen function, defined for E M as 

s'm{n9) 



(3.1) Cl2(^) = ^. 



n>l 



The Clausen function can also be defined as the imaginary part of Li2(e*^), 
where Li2(2;) is the dilogarithm function (see [21] for an extensive study of 
both dilogarithm and Clausen functions). 

Function Cl2(^) can be easily evaluated numerically. Definition (3.1) im- 
plies that for € M and n € Z 

Cl2(0 + 27rn) = Cl2(0), Ch{-e) = - Ch{9), C12{2tt - 6) = - ChiO) , 

and thus we only need to be able to compute this function for \9\ < vr, 
where it can be done very efficiently with the help of the following series 
representation (see [9] and formula 4.28 in [21]): 



Cl2(9)=3«-«ln( |«|( , 



where C{s) is the Riemann zeta function. Note that C(2^^) — 1 = 4~" + 9~" + 
• • • ; thus when |0| < vr the terms in the above series decrease as 0(n~^16~") 
and we have a very fast convergence rate. 

Theorem 2. Assume that a = ^ where m and n are coprime natural 
numbers. Define 

Q — f cot~^(cot(7rm/9) + (— l)™'"z™' sin(7rm/9)~^), if mp ^ Z, 
\ 0, ifmpGlj. 

Then for z > 

(j)(z) = exp( ^—(Chi29) - Ch(2Trmp) - Ch(29 - 27rmp)) 
V 2-Kmn 



(3.2) 



X (l + (-l)™"2cos(7rmp)z'"+z^'") 

n-l 



2m\-p/(2n) 



X 

fc=0 



n (1 + 2cOs(™(/5 + 2k + l))z" + ^2a)(n-2fe-l)/(2n) 
fc=0 

m~l , , . N (m-2i-l)/(2m) 

\[ + 2cos\^- {ap + 2j + I) y + z^j 



X 

i=o 
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The proof of Theorem 2 is based on the fohowing result: 

Lemma 1. If m and n are coprime natural numbers and |Im(2;)| < vr 

.m\ n z 
F\ z;i— 



n J m (-l)m«e"^ + 1 



n-l 



(3.3) +Y1 



Tri/n{n — 2k — 1) 



exp(z + TTim/n(2k + 1)) + 1 

k=0 

n ^-^--^ ■Ki/m{m — 2j — 1) 

m ^ — ^ 



m ^ exp(zn/m + iTin/m{2j + 1)) + 1 

Proof. The proof is in fact quite simple. The idea is to apply (2.8) 
and (2.9) in order to transform F{z;i^) into a sum of F{-;i), which can be 
evaluated explicitly. First we use (2.9) and obtain 

= - VFfl(z + 7ri(m-2i-l));-Y 
y n J m ^—^ \m n J 

Next, applying (2.8) to each function F{-,^) in the right-hand side of the 
above formula we find that 

Again, we apply (2.9) to each function F{-,in) in the right-hand side of the 
above formula and deduce 

(3.4) F z; i- = - V V F - + -(n - 2fc - 1) + -(m - 2i - 1); z . 

Now we only need to evaluate F(z;i), which can be done as follows: 

^, , e~^/^ f dx 
F{z;i) 



4 7]g cosh((3; — z)/2) cosh(a;/2) 
e~^/^ f dx 



2 cosh(x — z/2) + cosh(z/2) 

-z/2 ^ 



cosh(2;) -|- cosh(z/2) — 1 ' 

where in the last step we have used formula 3.514.1 from [17]. Finally, the 
double sum in the right-hand side of (3.4) can be reduced to single sums in 
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(3.3) using the following identity: 



1 1 ""^ 

1 + w"- I 



_|_ y,pTrim/n{ri—2k—l) ' 

fc=0 " 



which is just a different way of writing the partial fraction decomposition 
(2.14). □ 

Proof of Theorem 2. Using Proposition 1 and formula (3.3) we find 
that 

— ln((/)(z)) = (-1)™"+! ln(z™)mz™-i 



dz TTTJin 1 + (— l)™"'2cos(7rmp)2;™ + z^'"* 
p (-l)™"2cos(7rmp)mz™-i +2mz2™-i 

~ 2n 1 + (-l)'^'^2cos(7rmp)z'" + 

+ 2cos(7rQ(p + 2A; + l))az°'-'^ + 2az2a-i 

fc=0 
m— 1 

+ E 



2n / l + 2cos(7ra(p + 2A; + l))z" + 2;2a 
m-2j-l\ 2cos(7r/a(ap + 2j + l)) + 2z 



2m J l + 2cos{7r/a{ap + 2j + l))z + z'^' 



To complete the proof we just need to integrate the above identity and use 
the fact that (j){0) = 1. Most of the integrals are obvious, except for the 
following one: 

where G { — 1,1} and = cot~^(cot(/3) + ^2 sin(/3)~^). However, this integral 
is just a particular case of formula 8.18 in [21], and we refer the interested 
reader to this great book for all the details. □ 

Note that while Theorem 2 is different from the result obtained in [11], 
there are also some similarities. In particular, both results give explicit for- 
mulas for a dense uncountable set of parameters {a,p). If p is also a rational 
number of the form p = then 9 = 0, and formula (3.2) does not contain 
function CI2. It is easy to see that the set of points 

m k 
a= —,p= — 
n m 

is dense in the set of all admissible parameters A. It is probably true that 
in this case formula (3.2) can be reduced to expression (1.11) in [11]: due to 
the fact that the numbers m and n are coprime we can find integers k and 
I, such that k = In — km, thus p + k = l/a and the process Xt belongs to the 
class 



ON EXTREMA OF STABLE PROCESSES 



13 



4. Infinite product representation and analytic continuation. In this sec- 
tion we will study analytic properties of the Wiener-Hopf factor (j){z). The 
most important results are the representations in terms of the q-Pochhammer 
symbol and the double gamma function. 

Definition 4. For n G N we define the q-Pochhammer symbol as 

n-l 

{a;q)n = l[{l-aq'') 

and {a;q)o = 1. If \q\ < 1 we define (a;g)oo = nfe>o(l ~ o.q'')- 

Theorem 3. Assume that Im(Q!) > and define 
(4.1) g = e2™, g = e"2^*/". 

Then for \z\ < min{ y^M", \/|^} and |arg(2;)| < vr 

._{-zVqe-^'P;q)oo{-z''^e^'P'';q)oo 



(4.2) 



{-z^e^ip; q)oo{-z'^^e-^ip'^; q). 



Proof. Using Proposition 1 and series expansion (2.12) we find 



d 



5;lnW--)) = E 



A:>0 



1 



az 



a-l 



_j_ gTrjpqi— fc — 1/2 



■^a _j_ gvripag— fc— 1/2 
^a-1 



. _j_ g— Trip jj—fc— 1/2 



+ 



az 



_j_ g— TripoQ— fc— 1/2 



To complete the proof we only need to integrate both sides of the above 
identity and use the fact that = 1. □ 

As a corollary of Theorem 3 we can derive some results obtained in [11]. 

Corollary 1. If XtG Ck,i, then for |arg(z)| < vr 
r (z"(-l)l-'9(l-^)/2;9)fc 



(^(-1)1-^=^(1-0/2;^-), ' 

(^(-l)m-g(i+0/2.g),^, 
t (z"(-l)i+'g(i+^)/2;g)|,|^ 



ifl>0, 
ifl<0. 



Proof. Assume first that Im(a) > 0. We use (4.2) and the following 
identity for the q-Pochhammer symbol: 

(a; 9)00 



{aq^;q)c 



{a;q)n 
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(which can be easily obtained from Definition 4) to derive the above finite 
product representation. The last step is to let Im(a) — )■ 0"^ and use analytic 
continuation. □ 



A major disadvantage of Theorem 3 is that we can not extend it to the 
case that interests us, when a is a real number, since this would imply that 
l^l = kl = 1 the infinite products do not converge. Therefore our next 
goal is to transform (4.2) into an expression which is well defined for both 
real and complex values of a, and as we will see later, the double gamma 
function will be the magic key on our quest towards analytic continuation. 

Let us give some intuition on how one might guess that the double gamma 
function is in fact the right tool to use. Consider the following function: 



fH 



loo 



{-e°""^e~''^P°';q)c 



which enters into formula (4.2) with a change of variables z = e'^ . Using 
Definition 4 of the q-Pochhammer symbol we find that the numerator in 
the above expression for f{w) has simple zeros at points w = Wm,n + Trip for 
m > and n G Z, while the denominator has simple zeros aX, w = Wm,n + T^ip 
for m£7j and n < 0, where 

(4.3) Wm,n = Tri[{'^rn + l)a~^ + {2n+l)], m,neZ. 

We see that "half" of the zeros of the numerator coincide with "half" of 
the zeros of the denominator, and thus f{w) has simple zeros in the first 
quadrant of the lattice Wm,n + "tz/O (where m > and n > 0) and simple poles 
in the third quadrant Wm,n + Trip (m < and n < 0). Therefore it might be 
possible to separate the function f{w) into the ratio of two entire functions, 
if we could only find a function which has simple zeros in the given quadrant 
of the lattice. As it turns out, such a function was introduced by Alexeiewsky 
in 1889 and, luckily enough, all of its properties which are important for us 
were established by Barnes in [2] (see also [3]). 

We will follow the original paper [2], where Barnes introduces the dou- 
ble gamma function as an infinite product in Weierstrass's form: for all 
|arg(r)| < vr and z € C 

" ^az/T+bz^/(2T) 



G{z;t) = -e 

(4.4) 

nn'(i+ 



TJiT + n 

m>On>0 ^ 



X I I III 1-^ |g-2:/(mr+n)+2V(2(mr+n)2)_ 



Here the prime in the second product means that the term corresponding to 
m = n = is omitted. Note that by definition G{z; r) is an entire function 
in z and has simple zeros on the lattice mr + n, m<0, n<0. Barnes proves 
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that G{z;t) can also be expressed as a single infinite product of gamma 
functions 

(4.5) G{z;t) = ^e"^/^+^^'/(2-') TT r(mr) ,^^^^)^,y2^'imr) ^ 
tT(z) T(z + mT) 

^ ' ra>\ ^ ' 

where = S ^'^(''^('^)) digamma function (see [17]). The constants 

a and h are related to a and h as follows: 



d = a-jT, b = b+ ^ , 

where 7 = — ■0(1) is the Euler-Mascheroni constant. One of the most im- 
portant properties of the double gamma function is that it is quasi-periodic 
with periods 1 and r, that is, 

Giz + l-T) = T(^)G{z;T), 

(4.6) 

G{z + t;t) = (27r)("-i)/V-^+i/2p(^)(^(^. 
provided that constants a and b are chosen in a particular way 

(4.7) a = ^ln(2^r) + ^ln(r)-rC(r), 

(4.8) b = -Tln{T) -t^D{t). 

Here C(r) and D{t) are certain transcendental functions of r which can be 
computed as the following limits as m — t- -|-oo: 

3 

- ^^'(mr) + ^■0^^^ (rriT) + ©(m'^), 



D{t) = V'(^'^) + -tp'imT) ^^{mT] 

^—^ 2 T 

fc=i 

^3 



12^ ^ ^ 720 ^ ^ ^ ^ 
It turns out that with this choice of constants we also have G(1;t) = 1 
(see [2]). There exists a different and slightly simpler expression for these 
constants (see [20]), however we have decided to use the original Barnes 
formulas as they are more convenient for numerical calculations. It is also 
possible to give an integral representation for ln(G(z;r)) (see [20]) and sev- 
eral asymptotic expansions (see [6]). 
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The following two facts about the double gamma function will be very 
important for us later. The first result was derived in [2] and it is an analogue 
of the reflection formula for the gamma function 

(4.9) -27:iTG(l+z;T]G(l-z;-T^ - 



z 1 

T 



2 y v2 J {q;q)oo 

where q = e^'^*'^ and Im(r) > 0. The second result is the following transfor- 
mation formula: 

(4.10) G(z;t) = (2vr)^/2(i-i/r)^-.V(2r)+./2(i+i/r)-i^J^f . 

We were not able to find any reference for this result; however, it is not 
hard to prove it directly. It follows from the definition (4.4) that the zeros 
of G{z; t) coincide with the zeros of the function on the right-hand side of 
(4.10). Performing some straightforward but tedious computations one can 
check that the latter function satisfies both functional equations (4.6). Thus 
their ratio must be a function with two periods 1 and r, which is analytic 
in the entire complex plane, therefore assuming that r ^ M we conclude that 
this ratio must be a constant. The fact that the value of the constant is 
one can be established with the help of (4.6). We have included the detailed 
proof of formula (4.10) in the Appendix of the online version of this paper, 
see arXiv:1001.0991. 

Now we are ready to state and prove our main result in this section. 

Theorem 4. For a>0, p G (0, 1) and |arg(z)| < vr 

,G(l/2 + a/2(l + P + ln(z)/(7ri)); a) 



</>(z) = (2^V^)-"^ 



(4.11) 



G(l/2 + a/2(l -p + ln{z)/{TTi));a) 

G{l/2 + a/2(l + p- ln{z)/{Tri));a) 
G{l/2 + a/2{l- p-ln{z)/{TTi));a)' 



Proof. Assume that Im(a) > 0. We start with formula (4.2) and apply 
identity (4.9) to each q-Pochhammer symbol to obtain 

G(l/2 - l/(2a) - p/2 + ln(z)/(27ri); -1/a) 



G(l/2 - l/(2a) + p/2 + ln(z)/(27ri); -1/a) 



X 



G(1/2 + 1/(2q) +p/2 -ln(z)/(27ri) 



(4.12) 



G(l/2 + l/(2a) - p/2 - ln(z)/(27ri) 
G(l/2 + a/2 + pa/2 + aln(z) /(27ri) 



G(l/2 + a/2 - pa/2 + a\n{z) /(27ri) 
G(l/2 - a/2 - pa/2 - a\n{z) /{2m) 



G{l/2 - a/2 + pa/2 - a\n{z) /{2m] 



1/a) 



1/a) 
a) 



a) 
—a) 



-a] 
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Next, we transform the four double gamma functions in the first ratio (the 
ones having r = ±a~^) using identity (4.10) and after simplifying the re- 
sulting formula, we obtain (4.11). These calculations are straightforward but 
lengthy and tedious; the interested reader can find them in the Appendix of 
the online version of this paper, see arXiv:1001.0991. 

Now that we have established (4.11) for the case Im(a) > 0, we note 
that G{z]t) is analytic and well defined when r G M"*". Thus by analytic 
continuation (4.11) gives us an expression which is valid for a £ M"*". □ 

Now we can give a complete description of the analytic structure of (l){z). 



Corollary 2. Assume that a > and p G (0,1). The function 4>{z) 
can be analytically continued into the domain {z G C:|arg(2;)| <7r}, except 
that it has simple poles at points { — exp(±7ri(/9 — a~^)} when ap > 1. The 
function (j){z) has a branch point at z = while the function w i— )• (j){e^) is 
meromorphic in C and quasiperiodic with periods 2'Ki and 2TTia~^ 

1 _|_ paw+nia(l—p) 

(4.13) (/)(e"'+2"*) = ^ie'")—^—-jT^, 

^(g«;+2W") = 0(e"') ^ ., 

The function (j){e'^) has roots at 

{Wm,n + 7rip}m>0,n>0 U {Wm,n - '^ip}m<0,n<0 

and poles at 

{Wm,n - vrip}m,>0,n>0 U {Wm,n + '^ip}m<0,n<0, 

where Wm,n = 7ri[(2m + l)a~^ + (2n + 1)]. All roots and poles are simple if 
a is irrational. 



Proof. The analytic continuation result follows from expression (4.11), 
as the double gamma function G{z; r) has simple zeros at points z = — {mr + 
n), n,m > 0. The quasi-periodicity expressions and formulas for roots/poles 
of (^(e"') follow at once from the infinite product representation (4.2). □ 

Note that if {a,p) G A, then ap < 1, which implies that (j){z) is analytic 
in |arg(z)| < vr. This result coincides with analytic continuation result for 
(j){z) when Xt G Ckj, which was obtained in [11]. Also, Corollary 2 gives us 
an insight into the mysterious relation (1.1), which was the key to obtain 
explicit expressions for (l){z) in [11]. Assuming p = la~^ — k we find that all 
zeros and poles of 0(e"') lie on the same lattice Wm,n, and most of them 
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are canceled allowing us to express 4'{e'^) as a ratio of finite products of 
trigonometric functions. 

Applying the infinite product representation (4.5) to Theorem 4, we ob- 
tain the following expression for (j)(z) as an infinite product of gamma func- 
tions. The detailed proof can be found in the Appendix of the online version 
of this paper, see arXiv:1001.0991. 

Corollary 3. For a>0, pe {0, 1) and |arg(z)| < vr 

^ ^~ap/2^-ap{2C{a)+{a+l)D{a)) 



m>0 



1 a / \n(z) 
+ -{2m+l-p+ ^ ^ 



2 2V ' TTl 

„ , 1 a / ln(z] 
xr - + - 2m + l- p ^ ^ 



2 2V TTi 
r - + - 2m+l+p+ ^ ' 



2 2V ' vrf 



^ , 1 a f liviz] 
xr - + - 2m-hl + p ^ ' 



TTl 



-1 



5. Series representation and functional equations. As we have seen in 
the previous sections, rationality of certain parameters plays an important 
role in our results. In our next result there appear even more intriguing 
connections with Number Theory. It turns out that our next result is valid 
for almost all parameters a, except for rational numbers and for those real 
numbers which can be approximated by rational numbers unusually well. 

Definition 5. A real number x G M \ Q belongs to the set C if there 
exist b> 1 and c > such that the inequality 



P 

a 



c 



is satisfied for infinitely many coprime integers p and q. 

The set £ is a proper subset of the set of Lioville numbers, which satisfy 
the following, weaker property: for all p> there exists c> such that the 
inequality 



P 

X 



c 

< — 
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is satisfied for infinitely many coprime integers p and q. A celebrated result 
by Liouville states that any algebraic number is not a Liouville number, but 
this is also true for many other numbers. In fact, almost every number is not 
a Liouville number, as the set of Liouville numbers, while being dense in M, 
has zero Lebesgue measure (see Theorem 32 in [18]). The same is true for the 
set C: it is closed under multiplication by rational numbers, and therefore 
it is dense in R, yet it has zero Lebesgue measure. 

Theorem 5. Let a > and p e (0, 1). Assume that a ^ £ U Q. Then if 
\z\ < 1 and |arg(z)| < vr 



(p[z) = exp 

(5.1) 



sin(7rfcp) fc^fc 
^fcsin(7rfc/a) 



k>l 



^-^ fcsinfvr/ca) 

fc>i ^ ^ 



Proof. Assume first that z is a real number, such that \z\ < 1 and 
Im(a) > 0. Then Re(ib7ri/9 — ln(2)) > and formula (5.1) can be easily ob- 
tained using Proposition 1 and series expansion (2.13). The hard part is to 
extend validity of (5.1) to the case when a is real, and this is not a trivial 
matter. For example, when a is a rational number the series on the right- 
hand side of (5.1) is not defined, as some terms will include division by 
zero. 

Assume that Im(a) > and Re(Q;) is positive, irrational and does not 
belong to C Then using the above definition of the set C, we find that for 
every 6 > 1 there exists c> such that |/cRe(a) — n| > ckh~^ for all integers 
k and n. Therefore we have 

|sin(7rA;a)| > |sin(7rA;Re(a))| > ckb~^, k>l, 

which implies that the second series in (5.1) can be bounded from above by 

^-lg^Im(a)^^-l(^|^|Re(a))fc_ 

fc>l 

It is clear that the above series converges in the domain jzl^'^^") < b, |arg(z)| < 
vr. Now we take the limit Im(Q) — t- O"*" and use the Dominated Convergence 
theorem to conclude that (5.1) is true when z belongs to the above men- 
tioned domain and a is real, a ^ £ U Q. Since & > 1 is arbitrary we conclude 
that the series in (5.1) converges in the domain \z\ < 1, |arg(2:)| < vr. □ 

Theorem 5 was established independently in a recent paper by Graczyk 
and Jakubowski [15], where the authors applied series expansion to the in- 
tegrand in the formula for ■^\n{(p{z)), the latter was obtained from the 
Darling's integral. 
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We see that convergence of series (5.1) is intimately hnked with the degree 
to which we can approximate a number by rational numbers. The good news 
is that the series converges for |2:| < 1 for almost all values of a > 0, as the 
set £ has Lebesgue measure zero. However, we can easily exhibit a dense set 
of irrational numbers a for which the argument that we used to prove this 
theorem breaks down. For example, let us take an integer a > 1 and define 



— where qn+i = a''",qo = l. 



a = 

It is clear that if we take the sum of the first n terms we will have a rational 
approximation in the form Pn/Qm and the error of this approximation will 
be less then Ca~'^" for some C > 0. Then, assuming that p is not a Liouville 
number and considering a subseries of (5.1) corresponding to k = qm we 
obtain 



sin(7rg^ap) ^.m^a<,m > ^2 
qm sin(7rgma) ^ ^ 



(az") 



qr, 



m>l ^^^"^ 



and we see that this series cannot converge unless \az°'\ < 1, therefore the 
domain of convergence is strictly smaller then < 1. It is clear that we will 
have the same situation if we multiply a by any rational number, thus we 
have a dense set of real parameters a for which the domain of convergence of 
the series in (5.1) can be arbitrarily small. Of course the Lebesgue measure 
of the set of such "unlucky" values of a is zero, so Theorem 5 can still be 
used (with some caution) for numerical computations. 

When Im(a) > 0, one can obtain another series representation for (j){z) 
by applying the q-binomial theorem (see Theorem 10.2.1 in [1]) 

to (4.2). However, it seems to be very hard to say anything about the con- 
vergence of such series for a S M"*" and we did not pursue this further. 

In Corollary 2, we have seen that (j){e^) is quasi-periodic with two different 
periods. In the next theorem we collect other results on functional equations 
satisfied by 4>{z). 

Theorem 6. Assume that p G (0, 1). 
(i) For a> and |arg(2:)| < vr 

(5.2) ^(-^■,a,p] =z''f^{z;a,p). 
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(ii) Fora^{0,p ^) and |arg(2;)| < vrmmjl, ^} 
(5.3) (p{z;a,p) = 4>(z'^;-, 



,ap 



(iii) For n>2, a> n — 1 and |arg(2;)| < 7r(l — ^^^) 




(iv) For Q > and Im(z) ^ 



(5.5) 



z;a,l- p)4){-z] a, p) = {l + e 




where 5 = sign(Im(z)). 

Proof. Statement (i) follows from Theorem 4 and statement (ii) from 
the series representation (5.1). To prove (iii) we use Proposition 1 and The- 
orem 1, (ii), to find that for z S M 



To obtain (5.4) we apply transformation (5.3) in both parts of the above 
identity and rescale parameters. To prove statement (iv) we note that for 
{a,p) gA functions (j){—iz;a,l — p) and 4){iz\a,p) are the Wiener-Hopf 
factors <j)i{z) and </']~(z) thus functional equation (5.5) is just another way 
of writing the Wiener-Hopf factorization cj)'^ {z)(j)'^ {z) = (1 + '^{z))^'^. The 
result for general a > and p € (0, 1) follows by analytic continuation. □ 

Identity (5.2) first appeared in [14] and later in [15], while (5.3) is equiva- 
lent to Theorem 3 in [11] (it also appears in [14]). Note that both statements 
(i) and (ii) of Theorem 6 can be established directly with the help of Propo- 
sition 1 and results presented in Theorem 1, (i) and (iv). Another possible 
approach to establish these results is to use the same technique as in the 
proof of Theorem 3 in [11]: first we verify the transformation identity for the 
processes belonging to one of the Ck^i classes and then prove the general case 
by using the fact that classes Ck^i are dense in the family of stable processes. 
Note that Theorem 5 and functional equation (5.2) provide a convenient 
method to compute the Wiener-Hopf factors, as long as |z| 7^ 1 and a is not 
too close to being a rational number. 
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6. Mellin transform of the supremum functional. In the case of a general 
Levy process the Wiener-Hopf factor <i)q{z) gives us valuable information 
about the distribution of S'e(g), where e{q) ~ Exp((7), but the only way to 
translate this into information about the distribution of St is to perform a 
numerical Laplace or Fourier transform inversion in the g- variable. The scal- 
ing property of stable processes allows us to replace the integral transform in 
the g-variable with an integral transform in the z-variable, and a surprising 
fact is that the latter can be evaluated in closed form. The main tool in this 
section will be the Mellin transform of the supremum Si, defined as 

M{w)=M{w]a,p)=¥.[{SiT~^]. 

The function 7W(s) is well defined if Re(s) is sufficiently close to 1. Note 
that our M{s) corresponds to M~^{s) in [11]. Next, we define the Mellin 
transfom of (t){z) as 

/■oo r 

(t{s\a,p)= / z^~^4>{z]a,p)dz = / e^''4>{e^;a,p)dw. 



The link between M{s) and $(s) is the following identity, which goes back 
to [10] (see also [11]): 

(6.1) $(s;a,p) =r(s)r^l - ^M{1 - s;a, p). 

This identity can be easily established using the scaling property of stable 
processes, which implies that St = t^^'^Si 

poo 

^{s) = / z^-iE[e"^^'=(i)] dz = V{s)n{Se(i)V] = r(s)E[e(l)-^/"(5i)-^] 

JO 

= r(s)E[e(l)-^/°]E[(Si)-1 =r(.)r(l - ^)a^(1 - s;a,p), 

where we have also used the fact that e(l) ~ Exp(l) is independent of St- 

With the help of formula (6.1) we can translate functional equations for 
the Wiener-Hopf factor (i){z) into the following remarkable identities for the 
Mellin transform M{s). 

Theorem 7. The function Ai{s) can he analytically continued to a 
meroniorphic function in C A4{s) is quasiperiodic with periods 1 and a 

(0.2) M(s + i) = ^„„(,(,- i^))r(i - l)r(i - '-^)Mis), 

(6.3) M{s + a) = - sin(7r(ap - 1 + s))r(l - s)T{a - 1 + s)M{s), 

TT 
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and it satisfies 

A^/ N r(ap- 1 + s) r(l - p + (1 - s)/a) ...^ , 

r(i — s) r(i — (1 — s)/a) 
/ N ... N r(s) r(i + (i-s)/a) . / i-s 1 

i (z — s) i (1 — (1 — sj/aj \ a a 

Proof. First of all, we use (5.2) to find that (j){z) ~ z~^p as z — > +oo. 
Thus the function $(s; a, p) is well defined for < Re(s) < ap. Let us assume 
that a < ^ and define 

(6.6) G(s)= / —^-—e'"'^{e'")dw. 

Using the first functional equation in (4.13) we find the following identity: 
G(s)= I e"'^(/.(e"'+2-i)rfy;_e-2™p ^ e'"me')dw 

(6.7) 



/ e'""(/.(e"^+2'^')d?/;-e-2™''$(s). 



Corollary 2 and the condition a < | guarantee that ^{e™') is analytic for 
< Im(tt;) < Svr; thus we can shift the contour of integration M i— >■ M + 27ri 
in the integral in the right-hand side of the above identity and obtain 

e"'"0(e'"+2^^) du; = I (\>{e'")dw = e-'^''^'^{s). 

Jm. 

Combining (6.7) with the above equation we conclude that 
(6.8) G{s) = (e-2"^ - e-2™'')$(s). 

The integral which defines function G{s) converges for < Re(s) < a{l + p), 
and therefore (6.8) allows us to continue the function $(s) analytically into 
a larger domain < Re(s) < a(l + p). Next, one can use definitions of $(s) 
and G{s) to check that 

G{s) + e™(^+'')G(s + a) = (1 - e-^^^^O^C^), 

and combining the above two equations we obtain 

^/ N sin(7rs) ^, . 

$ s + a) = - . , , ^ T\T^ s • 

sm(7r(s + a(l — p))) 

Formula (6.3) follows easily by applying (6.1) to the above identity. Equation 
(6.8) extends ^>(s) analytically into the domain < Re(s) < a(l + p), which 
means that Al(s) is well defined for 1 — a(l + p) < Re(i(j) < 1 and applying 
(6.3) repeatedly we can extend A^(s) to a meromorphic function in the 
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entire complex plane. The restriction a < 5 can be removed by analytic 
continuation to other values of a and the second identity (6.2) can be proved 
similarly by starting with the second functional equation in (4.13). 

To prove (6.4) and (6.5), we apply the Mellin transform to both sides of 
equations (5.2) and (5.3) to obtain 

^>(s;a,p) = ^>(a/)- s;a,p), ^>(s;a,p) = -$( -;-,ap ). 

a \OL a J 

Functional equations (6.4) and (6.5) follow from the above identities and 
(6.1), the details are left to the reader. □ 

A surprising fact is that quasi-periodicity of M{s) allows us to find an 
expicit formula for this function, and having already developed all the nec- 
essary tools we can enjoy a rather simple and straightforward proof. 



Theorem 8. For s G 
M{s) = a 



s-i G{ap;a) 



Gal-p+l;a 
(6.9) 

G{a{l-p) + 2-s;a) G{a-l + s;a) 

G{ap — l + s;a) G{a + l — s;a) 

If Xf £ Ck,i and I > 0, then 

M(<i) r(g) ^ sin(7r/a(g - 1 + j)) 

r(l-(l-s)/a) y sin(7rj7Q) 

(6.10) 

k 

nsm(7rai) 
j 

If Xt £ Ck,i and I < 0, then 



^ sin(7r(l — s + aj)) 



\k\ — l 

j^l^^^ ^ r(l + (l-s)/a) -j-r sin(7r(g- l + «j)) 
r(2 — s) sin(7raj) 

(6.11) 

^ TT smjirj/a) 

sin(7r/a(l — s + j)) 



Proof. Formula (6.9) can be established using the classical approach 
from the theory of elliptic functions. Let us denote the function on the 
right-hand side of this formula as H{s). Using quasi-periodicity properties 
of G{z;t) [see formulas (4.6)] we find that II{s) satisfies (6.2) and (6.3) 
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(the detailed computations can be found in the Appendix of the onhne ver- 
sion of this paper, see arXiv:1001.0991). This imphes that Ai{s)/ H{s) is a 
meromorphic function which is periodic with periods 1 and a. An integer 
hnear combination of periods is also a period; thus, assuming that a G M"*" is 
irrational we can find arbitrarily small periods of the form ma + n, and there- 
fore the function M{s)/H{s) must be constant. The value of this constant is 
equal to one since A4(l) = H{1) = 1. Formulas (6.10) and (6.11) can be ver- 
ified using exactly the same approach, again the detailed proof can be found 
in the Appendix of the online version of this paper, see arXiv:1001.0991. 
Otherwise, formulas (6.10) and (6.11) can be established directly, using the 
fact that ap = I — ka and transforming (6.9) with the help of the functional 
equations for the double gamma function (4.6). Finally, we extend the result 
to the case when a is rational by analytic continuation (note that all the 
formulas in Theorem 8 are well defined when a G Q). □ 

The Mellin transform of the supremum was evaluated for the three special 
cases: Co,i (spectrally negative process) and Ci^i, Ci_2 (processes with two- 
sided jumps) in Corollary 2 [11], but they do not agree with our results 
presented in Theorem 8. It seems that formulas (3.4), (3.5) and (3.6) in [11] 
are not correct, as they violate identity = 1 (which follows from the 

definition of the Mellin transform). 

It is interesting to note that formula (6.9) is valid even in the limiting 
case a = 2 and p = 1/2, that is, Xt is Brownian motion \plWt- Using the 
properties of the double gamma function given in (4.6) one can check that 
(6.9) reduces to 

Mis) -2^-^^-^ 

^^'^-^ r(i/2)' 

and the above formula can also be verified directly using the distribution of 
the supremum of Brownian motion (see [8]). 

The results in Theorems 7 and 8 are surprisingly explicit, and later in 
Section 7 we will use them to derive series representations and asymptotic 
expansions for the probability density function of S\. Here we also present 
two immediate corollaries of Theorems 7 and 8. 

Corollary 4. Assume that a G (^, 1) U (1,2). Let Xt and Xt be stable 
processes with parameters {a,p) and {^,ap) having corresponding supremum 
processes St and St- Then we have the following identity in distribution: 



(6.12) £1 



Si 

^2 



5i 

£4 



where £i ~Exp(l) and all random variables are assumed to be independent. 
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Proof. To prove this statement we change the variable in (6.5) s i— t- 
1 — a(l — s), rewrite the resulting identity as: 



r(l + a - as)T{s)M{l -a + as; a, p) = T{l-a + as)r{2 - s)M 



and use the following facts: (i) r(s) is the Mellin transform of the Exp(l) 
random variable; (ii) if M{s) = K[X'^~^] is the Mellin transform of the ran- 
dom variable X, then M(l — o + as) is the Mellin transform of y = X'^; (iii) 
the Mellin transform of the product of two independent random variables is 
equal to the product of their Mellin transforms. □ 

Corollary 4 should be compared with Theorems 12 and 13 in [26], which 
describe identities in distribution satisfied by products of powers of stable 
random variables, and with Theorems 1.1, 1.2 in [14]. Also, note that using 
the scaling property of stable processes we can rewrite identity (6.12) in a 
more symmetric form 

5*6(1) " _d Se{l) 

ei £2 



Corollary 5. If Xt & Ck,i and I > 0, then Si satisfies the following 
identity in distribution: 

1/q 



.13) 



5i X 



i-i 



7{a-ii} 



■-£2 



7l-{ai} 



where {x} G [0, 1) denotes the fractional part of x, Si ~ Exp(l), 7^ denotes a 
gamma random variable with E[7fc] = Var[7fc] = k and all random variables 
are assumed to be independent. 

Proof. The proof is very similar to the proof of Corollary 4: we rewrite 
(6.10) as 



i-i 



sm{7rj/a) 



L sin(7r(l — s + aj)) ' 



M{s)t(i-^-^]Y\- . , , , ^ 

V " J j-J^sm{TT/a{s - 1+ j)) ^_ 

and use the fact that for < a < 1 the Mellin transform of y = 7a/7i- 
(where 7^ and 7i_a are independent) can be computed as 

r(a + s-l)r(2-a-s) 



E[y^"i] = E[(7a)^-']E[(7i^a)'' 
sin(7ra) 
sin(7r(a + s — 1)) 



r(a) 



r(i-a) 



□ 
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Corollary 5 should be considered as a generalization of the corresponding 
result for Brownian motion, that is, when Xt = ^/2Wt. In this case we have 

o 1/2 d 
Si X =£2. 

This fact can be established by applying the scaling property to the identity 
S'e(i) ~ Exp(l), which follows from the Wiener-Hopf factorization for Brow- 
nian motion. Identity (6.13) also can be rewritten in the form Sti = T2 with 
an obvious choice of random variables Tj. It might be possible to establish 
a similar identity in the general case, using formula (6.9) instead of (6.10). 
However, it is not clear how to connect the double gamma function with the 
Mellin transform of some random variable and we will leave this for future 
work. 



7. Probability density function of the supremum functional. In the last 
two years there have appeared several interesting and important results re- 
lated to the density of the supremum = supfX^ : < s < 1}. In the spec- 
trally positive case Doney [12] has obtained the first asymptotic term for 
the density 

p(x) = -^F(Si < x) 
ax 

as X — )• +00; Bernyk, Dalang and Peskir [4] have derived an explicit conver- 
gent series representation for p{x) and Patie [23] has obtained a complete 
asymptotic expansion of p{x) as x — ?• +00. In the case of a general stable 
process Doney and Savov [13] obtain the first term of asymptotic expansion 
of p{x) as X — )• 0"*" or X — )• +00, and they also mention that it is possible 
to obtain higher order asymptotic terms as x — t- +00 (see Section 5 in their 
paper). See also [25] for asymptotic results on distributions of functionals 
of a certain random walk related to stable processes and [24] for an explicit 
infinite series representation for the density of the first hitting time of a 
point in the spectrally positive case. 

The following two theorems summarize our main results in this section. 

Theorem 9. Assume that a ^ Q. Define sequences {am,n}m>o,n>Q o.nd 

{bm,n}m>0,n>l 0,S 



a. 



m,n 



(7.1) 



r(l — p — n — m/a)T{ap + m + an) 

sm{7r / a{ap + j — 1)) yj sin{TTa{p + j — 1)) 



m 

■j--^ sin(7r j'/a) -j--^ sin(7rQj) 



, , T{1 — p — n — m/a)T{ap + m + an) 

(7.2) Om,n = WT-, ; ; TTFT ^ 
1 (1 + n + m/a)i. (— m — an) 



m,n- 
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Then we have the following asymptotic expansions: 

(7.3) J] j;a^,„x'"+"", x^0+, 



n>Om.>0 



(7.4) pix) ~ Yl E bm,n+lX-", X ^ +00. 

n>Om>0 

Remark 1. Asymptotic expansion (7.3) (and other similar expressions) 
should be understood in the following sense: for any c > we have 

p(x) = x"^-i a™,„x'"+'^" + 0(a;"^-^+'=), x^0+. 

m,n>0 
m+an<c 

Theorem 10. Assume that Xt G Ck,i. If I > 0, then for n G {0, 1, . . . , fc} 
and niGZ, we define 

_^ (^_l-^m{k+l)+nl+l 

^"''^ ~ r(l + n + m/a)T(-m - an) 

(7.5) 

-i-j- sin(7r/Q!(j + m)) -r-j sin(7ra(j + n)) 
■j^ sin(7rj/a) sin(7r aj) ' 

while if I <0, then for m € {0, 1, . . . , |/|} and n G we define 

(^_l-^mk+n(l+l)+l 

^""'^ ~ r(l + n + m/a)r(-?n - an) 

(7.6) 

A: —1 / — m 

T-r sin(7rQ;(j + n)) -pr sin(7r/a(j + ?n)) 
-'^ sin(7rQ:j) -j^ sin(7rj'/a) 

T/ien if a £ (0, 1) and / > {/ < 0} we /ia?;e a convergent series representa- 
tion 

k 

n=l m>0 

(7.7) 



s p(x) = - E E ^"1." 

t m=On>l 



and an asymptotic expansion 

k 



n=Om<—l 
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(7.8) 



I m=ln<k 



Similarly, if a £ (1,2) and l> {I < 0}, we have a convergent series repre- 
sentation 



k 

m—an—1 



(7.9) 



n=Om<—l 



(. m=ln<k 



and an asymptotic expansion 

k 



n=l m>0 

(7.10) 

1^1 



r w 

I p{x) ~ - X] X] 

I m=On>l 



m-an-l\ x ^ +00. 



Infinite series in (7.7) and (7.9) converge uniformly on compact subsets of 
R+. 

First we need to establish some technical results which describe the ana- 
lytic structure of 7W(s). 

Lemma 2. The function Ai{s) can be analytically continued to a mero- 
morphic function in C If ct and Xt ^ C^^i for all k and I, then M{s) 
has simple poles at 

{^m,,n}m>l,n>l = {m + an}m>l,n>l, 
{^m,n}m>0,n>0 = {I — ap — m — an}m>0,n>0, 

with residues 

(7.11) Res(A^(s) :s+ „) = -6m-i,n, Res(A^(s) : s;; „) = a^^^. 

In the case when Xt G Ck^i and I > {I <0}, the function Ai{s) has simple 
poles at Sm,n = m + an, where 

m<l-l, n G {0, 1,2, . . . , A;} or m > 1, n G {1, 2, . . . , A;}, 

{m E {1,2, 3, ... , |/| + l},n > 1 or m e {2,3, ... ,\l\ + l},n< k}, 
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with residues 

(7.12) Res(7W(s):s„,„) = c+_i „, {Res(A^(s) : = c-_i_„}. 

Proof. Iterating identities (6.2) and (6.3), we find 
M{1 + s + m + an) 
= M{l + s)a"'+''' 

^ r(l - (1 + g + j)/a)r(l + is + j)/a) 
- - {s + j)/a)T{p +{s + j)/a) 

'i-r^ T[—s — m — aj)T{a + s + m + aj) 
^_}-^ r(l — ap — s — m — aj)T{ap + s + m + aj) 

The right-hand side of the above equation has a simple pole at s = 0, which 
comes from r(— s — m) (take j = in the last product). The residue of 
r(— s — m) at s = is equal to (— l)™~^/m!, and therefore 

Res(A^ (s) : 1 + m + an) 

(-l)'^-ia'^+"r(a + m) 



(7.13) 



m!r(l — ap — m)T(ap + m) 
^^YT r{l-{l+j)/a)r{l+j/a) 

f}o r(i-p-jya)r(p + j/a) 

n — 1 

Yj r{—m — aj)T{a + m + aj) 
r(l — ap — m — aj)r{ap + m + aj) 

The above expression is equal to —bm^n, which can be verified with the help 
of the reflection formula for the gamma function. The detailed computa- 
tions can be found in the Appendix of the online version of this paper, see 
arXiv:1001.0991. The value of the residue at s~ „ can be easily established 
with the help of the reflection formula (6.4). Finally, when Xt G Cfc,i equation 
(7.12) can be derived by computing the residues of Ai{s) given by formula 
(6.10). Otherwise one could establish it as a corollary of the previous result 
by checking that am,n = c^i_m,fc-n ^^'^ ^^.^ = -Cm,n depending on ±1 > 0. 
The latter approach with detailed calculations can be found in the Appendix 
of the online version of this paper, see arXiv:1001.0991. □ 

Lemma 3. For x gM. we have as y ^ oo, y G M 

hi{\M(i: + iy)\) = -^(a(l - p) + 1 - ap) 

(7.14) ^; 

+ o{y). 
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Proof. Equation (4.5) in [6] gives us the following asymptotic expan- 
sion valid as z — 7- oo in the domain |arg(z)| < vr: 

HGiz;a)) = l^Hni^) - - ^.ln(.) 

(7.15) +l(l±M5!)+in(2vra) + l)z 

where c(a) is some constant depending on a. The asymptotic formula (7.14) 
can be obtained from (6.9) and (7.15) using the following asymptotic expan- 
sion for the logarithm function: 

A 

ln(^ + s) = ln(s) H ir^ + 0{s~^), s-^oo, |arg(z)|<7r, 

s zs^ 

after some straightforward but tedious computations, which can be consid- 
erably simplified with the help of symbolic computation software. □ 



Proof of Theorem 9. Equation (7.14) and the fact that ap<l for 
all (a,/o) G A imply that A4(s) decreases exponentially as Im(s) — >• oo; thus 
5*1 has a smooth density function p{x), which can be obtained as the inverse 
Mellin transform, 



(7.16) 



p{x) = — [ M{s)x ""ds. 



By shifting the contour of integration 1 + iM i— )• c + iM where c < and taking 
care of the residues at points s~ „ we obtain 



1 



X X 



H / M{s)x ""ds, 

2vr^ Jc+iR 



(7.17) p(x) = J^Res(A^(s):s^,J 

where the summation is over all m > 0, n > 0, such that s^^> c. The inte- 
gral in the right-hand side of (7.17) can be estimated as follows: 



— [ Mis)x-'ds 
2m 7c+iR 



< 



x 

'2^ 

x"^ 
~2^ 



Mic + it)x~'^ dt 
\M{c + it)\dt = 0{x^ 



This establishes the asymptotic expansion (7.3). The proof of (7.4) is iden- 
tical, except that now we need to shift the contour of integration in (7.16) 
in the opposite direction. □ 
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Proof of Theorem 10. The asymptotic expansions (7.8) and (7.10) 
can be derived in the same way as in the proof of Theorem 9 (or obtained as 
its corollaries); thus we only need to establish convergence of (7.7) and (7.9). 
Let us assume that a G (1,2) and I > 0. Choose cq G (0, 1) such that cq ^ 
m + an for n G {0, 1,2, ... ,k} and m G Z. We take to be a large positive 
number and shift the contour of integration in (7.16) 1 + iM i— t- cq — + 
while taking into account residues at Sm,m hence 

(7.18) p(x) = ^Res(X(s):Sm,„,) X X"' 



1 



2iri 



where the summation is over m,n such that cq — N < Iie{sm,n) < 1- 

Using (6.10) and the reflection formula for the gamma function we find 
that for some constant C G M 



(7.19) 



Mis) = C 



r((l - s)/a) n5=oSm(vr/a(s - 1 + j)) 



r(i 



n)=o sin(7r(l-s + aj)) 



Now we need to prove that as — t- +oo the integral in the right-hand side of 
(7.18) converges to zero for all x G M''". Intuitively this is clear, since the ratio 
of gamma functions in (7.19) decreases to zero faster than any exponential 
function as Re(s) — t- — oo, while the other factor is just a ratio of periodic 
functions in s. The rigorous proof can be obtained as follows: 



/ 



M{s)x "ds 



Ico-N+iR 

where we have denoted 



<Cx 



T{iN -co + it)/a) 



T{N-co + it) 



g{t)dt. 



g(t) = e'^/"'!*! ]^|cosech(7r(t + i{aj - co)))|. 

j=0 

Using Stirling's approximation for the gamma function one can check that 
for all a; > 0, the function 



X 



^_^r{{N -co + it)/a) 
r{N-co + it) 



converges to zero as — )• +oo (uniformly in t G M); thus the integral in the 
right-hand side of (7.18) vanishes as A^ — >• +oo, and we have a convergent 
series representation (7.9). The convergence of series (7.7) can be established 
in the same way, except that now we have to shift the contour of integration 
in the opposite direction. □ 



It is important to note that all the asymptotic expansions and series rep- 
resentations for p{x) presented in Theorems 9 and 10 can be differentiated 
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N times term-by-term, where > 1 is an arbitrary integer. For the series 
representation this foUows easily by the standard argument of interchang- 
ing derivative and summation (both the series and its derivatives converge 
uniformly on compact subsets of More work is needed to establish 

a similar result for the asymptotic expansions, as it is not generally true 
that one can differentiate asymptotic expansions term-by-term. A classic 
counter-example is provided by the function 

f(x) = — \- cos(e'^^) x~"- asx— >+oo, 

1 — X ^-^ 

n>0 

for which the asymtptotic expansion of f'{x) cannot be obtained by simply 
taking term-by-term derivative of the asymptotic expansion of f{x). Thus 
in order to establish the result on the asymptotic expansion of p(^) (x) we 
would have to repeat the steps of the proof of Theorem 9: first we take 
A^th derivative of both sides of (7.16). Then we interchange the order of 
integration and differentiation in the right-hand side of (7.16); this can be 
easily justified using the standard uniform-convergence argument. Finally 
we shift the contour of integration to obtain asymptotic estimates. 

It is very likely that the asymptotic expansions given in Theorem 9 can 
also be interpreted as convergent series; however, it seems to be very hard 
to prove this fact analytically. The above proof of Theorem 10 was based on 
two facts: M{s) decays faster than any exponential function as Re(,s) — )• ±00 
[the sign depends on whether a G (0,1) or a G (1,2)], and that M^{s) is 
essentially a product of a function which decays very fast and a function 
which is periodic. The first fact is still true in the general case: using (7.15) 
and (6.9) one can prove that as Re(s) — >• 00 and < e < |arg(s)| < vr — e, we 
have 

log{M{s)) = {s-l) Hs) (1 - ^) + 0{s), 

log{M{-s)) = -{s + l) ln(s) (1 - ^) + 0{s). 

The major problem now is that when Xt ^ Ck,i we do not have any peri- 
odicity. Moreover, to make matters worse, when we move the contour of 
integration farther away from zero, the poles of Ai{s) become more and 
more dense, and it is very hard to find an upper bound on |7W(s)| for small 
values of Im(s). 

One can also see that the behavior of the coefficients CLm.n and bm n is much 
more unpredictable compared to Cm,n, and it is hard to say anything about 
the growth/decay of these coefficients as m or n becomes large. Numerical 
results, however, indicate that the product of ratios of sin(-) functions in 
(7.1) remains bounded as m or n becomes large, as long as a is not too close 
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to a rational number. Thus it seems reasonable to expect that the following 
conjecture is true: 

Conjecture: Assume that a ^ £ U Q. If a G (1,2) {a £ (0, 1)}, then the 
infinite series (7.3) {(7.4)} converges to p{x) for all x > 0. 

Acknowledgments. The author would like to thank two anonymous ref- 
erees for many detailed comments and constructive suggestions. 

SUPPLEMENTARY MATERIAL 

Appendix A: Detailed proofs of some results related to the double gamma 
function (DOI: 10.1214/10-AOP577SUPP; .pdf). This supplement material 
provides detailed computations needed to derive formulas (4.10), (4.11), 
(7.1), (7.2), (7.5) and to prove Corollary 3 and Theorem 8. 
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